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Abstract 

The zero range process is of particular importance as a generic 
model for domain wall dynamics of one-dimensional systems far from 
equilibrium. We study this process in one dimension with rates which 
induce an effective attraction between particles. We rigorously prove 
that for the stationary probability measure there is a background 
phase at some critical density and for large system size essentially all 
excess particles accumulate at a single^ randomly located site. Us- 
ing random walk arguments supported by Monte Carlo simulations, 
we also study the dynamics of the clustering process with particular 
attention to the difference between symmetric and asymmetric jump 
rates. For the late stage of the clustering we derive an effective master 
equation, governing the occupation number at clustering sites. 

Key words: zero range process; nonequilibrium phase transition; 
equivalence of ensembles; relative entropy. 

1 Introduction 

Low dimensional stochastic particle systems far from equilibrium have a 
much richer structure than their equilibrium counterparts. In particular, 
even in one spatial dimension there is the possibility of a phase transition. 
On the other hand, we do not have available general criteria which would 
allow us to determine, for example, the phase diagram. A recent attempt in 
this direction is a proposal by Kafri et al. fT], who study phase separation 
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in one dimension. Roughly speaking, they map the domain wall dynamics 
of clusters to a zero range process, for which analytical tools are available. 
Under suitable conditions, the clusters tend to grow. This coarsening is an 
intriguing phenomenon already on the level of the zero range process itself, 
and we will investigate it in much greater detail than available so far. 

The zero range process is a stochastic particle system on the lattice Z*^ 
where the jump rate g{k) of a given particle depends only on the occupation 
number k at its current position. This model was originally introduced as 
a simple example of an interacting Markov process |2]. Various properties 
have been established, among them the existence of the dynamics under very 
general conditions, classification of invariant measures, and hydrodynamic 
limits [3 in E]- If g{k) is decreasing in k, then this induces an effective 
attraction between particles, as first noted in jniEj. As a result, there is a 
critical background density and excess particles condense on a non-extensive 
fraction of the volume. 

Such a clustering phenomenon can be studied on two distinct levels. 
Firstly, the phenomenon is present already in the steady state. Based on 
some results for large deviations of independent, identically distributed ran- 
dom variables [HI IH| and under general assumptions on g{k), we will prove 
that for a typical steady state configuration there is a background phase 
at some critical density. Any additional mass is concentrated at a single, 
randomly located site. For the model introduced in [2,, we will analyze the 
statistical properties of the background phase in detail. Secondly, there is 
the dynamics of clustering with the steady state appearing only in the long 
time limit. In a very recent preprint Godreche ^U] addresses this problem. 
He assumes a uniform initial density and investigates numerically how the 
probability distribution of the number of particles at some given site evolves 
in time, with particular attention to the macroscopic component of that 
distribution. In view of our static result we study the dynamics of Evans' 
model [7] , pursuing a somewhat different approach: During the initial nucle- 
ation process random sites are selected, at which a macroscopic number of 
particles accumulates. We investigate the effective dynamics of the number 
of particles at such cluster sites, in particular how the smaller occupation 
numbers become extinguished to the benefit of the larger ones. 

Our contribution should be understood study in the context 

of phase transitions in one-dimensional systems far from equilibrium, which 
has been a topic of major interest in the past decade (see [TT| IT^ and refer- 
ences therein). Of particular interest is the occurrence of phase separation in 
systems with two conservation laws (see e.g. ^3) and jlllllS]) whose macro- 
scopic behavior has only recently been examined to some extent fOl El E| • 
Given the correspondence to the zero-range process our results also provide 
new information on the stability of domain walls (shocks) which separate 
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macroscopic regions of different phases in two-component systems. Domain 
wall stability already proved to be a key ingredient in the theory of boundary- 
induced phase transitions in systems with one conservation law |19| 1201 1^ 
and thus may shed light on boundary-induced spontaneous symmetry break- 
ing 122] in two-component systems. 

To give a brief outline: in the following Section we discuss the stationary 
measures for the zero range process and in Section 3 the two main results on 
the equivalence of ensembles and the structure of the condensed phase are 
established. In Sections 4 and 5 we study the condensation transition for 
the model introduced in first on the level of the stationary measure, and 
secondly through the dynamics of the clustering. 

2 Zero range process and its infinite volume 
stationary measures 

For notational simplicity we restrict ourselves to a description in one space 
dimension, but our results on the equivalence of ensembles hold for arbitrary 
dimension. Rather than defining the zero range process directly on an infinite 
lattice (cf. El), we first consider a finite system, compute the (unique) 
stationary measure and analyze it in the limit of infinite system size. 

Let us consider a zero range process on the one-dimensional lattice = 
{!,..., L} of L sites with periodic boundary conditions. Let r^^. G N be 
the number of particles on site x G A^,, where N = {0, 1, 2, . . .}. The state 
space is given by fli = and we denote a particle configuration by rj = 
{Vx)x(^Al ^ ^L- At a given site x G A^, the number rj^ of particles decreases 
by one after an exponential waiting time with rate g{rix) and the leaving 
particle jumps to site x + y with probability p{y). The jump probabilities 
p : Z ^ [0, 1] are normalized, '^yPiy) = 1, p(0) = 0, assumed to be of finite 
range, pijj) = for \y\ > R, and irreducible, p(l) > 0. For the dynamics 
to be well defined in the limit L — > cxd and to be nondegenerate the rate 
function : N — > [0, oo) has to satisfy 

sup|^(A; + l) -^(A;)| < oo , ^(A;) > ^(0) = for all A; > 0. (2.1) 

km 

The generator of the zero range process is then given by 

L R 

i^fM = E E 9{v.)piy){firi'''''^') - firi)), (2.2) 

x=l y=—R 

regarded as a linear operator on C(i7L,]R). We used the shorthand rj^^^^y = 
Tjz — 5{z,x) + 5(^z,{{x + y — l)modL) + l) for all z G A^, where 5(., .) denotes 
the Kronecker delta function. The process conserves the number of particles 
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Sl(^) = J2t=i''lx, thereby dividing the configuration space into the finite, 
invariant subsets ^l,n = {v ^ ^Ll^Liv) = ^} with N ^N. Restricted to 
C{Ql,n, K) with L, N fixed, £ is a finite dimensional matrix and the process 
is well defined. However for L = oo this is true only for "reasonable" initial 
conditions and under the assumption (j2.1|) . see jHlin]- 

The following results for the stationary measures are well known and 
taken from [21 El El- The zero range process ()2.2|) on Qn,l is an irreducible 
Markov jump process with the unique stationary measure 



(2.3) 



The weight W is given by 



H--W:=n^ (2.4) 
and the normalizing partition function is 

L 

Z{N, L)= W{r^,) 5(S^(r/), N) . (2.5) 

77eQi x=i 

Clearly ()2.3p resembles a canonical ensemble in statistical mechanics. There- 
fore, in the limit of large system size L, N ^ oo with fixed particle density 
p = N/L, ()2.3j) is expected to be equivalent to a grand canonical product 
measure, which is defined through 



x=l 



and where the fugacity > is adjusted to fix the average density. 

Let (pc be the radius of convergence of the grand canonical (one site) 
partition function 

CO 

Z{<p) = Y,W{k)<p'. (2.7) 

k=0 

The measure ()2.6|) is well defined for fugacities G [0, (pc) and its average 
particle density p(0) function of is given by 

p(^)=f:tM*)=^^^. (2.8) 

A:=0 ^ 

The range of p is the interval [0,pc), with p(0) = and pc = lim^yif,^ p{(j)) 
the critical density, (p p{(j)) is strictly increasing and we denote the inverse 
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function on [0, Pc) by 0(p). If 0c = oo, then = oo (see jl], Lemma 2.3.3), 
whereas for 0^ < cc, both p^ = oo and pc < oo are possible. In the second 
case Z{(j)c) < oo (see j;4j, Lemma 2.3.3) and z/^^ is a well defined probability 
measure with {rjx)v^^ = Pc- Thus we set 

r inverse of p(0) , for p < p 
\ 4>c , for p > Pc ' 

where pc can be either finite or infinite, cf. Figure 1 in Section 4. The reason 
for this particular convention will become clear in the next section. In this 
way we may regard the measure z/^ also as a function of p through 

The link between canonical and grand canonical measures is given through 
the pointwise limit of the n-point marginal, 

n 

lim/.lf^l'^(k) = n^0(p)(fc.)- (2.10) 

i=l 

Here [a] denotes the integer part of a G M and for every n G N, x = 
(xi,...,x„) G with Xi 7^ Xj for i j, k = {ki, . . . , kn) G N" the n- 
point marginal is defined as p^'^(k) := fJ^^'^{{r]xi = ki,...,rix„ = kn}). 
Since the measure p^'^ is permutation invariant (see ()2.3|) ). the marginals 
do not depend on the sites x individually, but only on their number n. For 
this reason only n is specified in our notation. 

A rigorous result on this equivalence is available in Appendix 2, but 
it does not cover the supercritical case pc < oo and N/L > pc. It will be 
discussed in Section 3, where particular attention is given to the statistical 
properties of the excess density. If pc < oo, then z/^^ is well defined and by 
necessity decays subexponentially, which will be explained in (|3.12|) . Exam- 
ples with power-law decay are given in jB], [Ij (Example 2.3.4). In [7j Evans 
studies the dependence on the power-law exponent in more detail and intro- 
duces a generic model where pc can be either finite or infinite, depending on 
a system parameter. Its stationary and dynamical properties will be studied 
in detail in Sections 4 and 5. 

The jump probabilities p do not influence the stationary measures, but 
they play an important role for the relaxation dynamics, see Section 5. The 
zero range process is reversible if and only if p is symmetric. The stationary 
current is given by j = ra{p){g{r].j)) ■, where ni{p) = J2y=~RyPiy) 
first moment of p, which for non-symmetric jump probabilities is generically 
non-zero. The average jump rate (g) of the two ensembles is given by 

Z{N-1,L) . 
{9{Vx))t,^,L = (canonical), 

oo 

{9iVx))u^ = ^9{k)v<i,{k) = (grand-canonical). (2.11) 

k=0 
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Thus for m{p) > the grand canonical current j(p) = m{p)(f){p) is monotone 
increasing in p, approaching its maximum value m{p)(j)c as p ^ pc, and 
correspondingly for m{p) < 0. 



3 Equivalence of supercritical measures 

We consider the zero range process with p^ < oo and supercritical canonical 
measures with p = N/L > p^. The heuristic picture, developed in [6^ and 
is that most sites of the system are distributed according to i/<^^ with mean 
occupation number pc. For large L, the {p — pc)L excess particles presumably 
condense on a few sites. If so, locally one will observe the grand-canonical 
ensemble with (j) = (pc- This picture is made precise in 

Theorem 1. (Equivalence of ensembles) 
Let p (f){p) be defined as in ()2.8|) . ()2.9|) . Then for every p G [0, oo) the 
n-th marginal has the pointwise limit 

n 

lim/.|f^l'^(k) = n^^<^(p)(^.)- (3.1) 

i=l 

The canonical partition functions converge as 

hm j\ogZ{[pL],L) = logZ(0(p)) -plog0(p) . (3.2) 



Proof. It is convenient to characterize the distance between the canoni- 
cal and grand-canonical measures through the relative entropy S. For two 
arbitrary probability measures p, on a countable set f2 it is defined as 

r V p(cu)log44 , if p < 
S{p\v) = \^n ' "^("^ ^ (3.3) 

[ oo , otherwise 

where p <^ z/ means that p is absolutely continuous w.r.t. v. It has the 
properties (see e.g. [2^) 

5(p|z/)>0, S{p\u) = Q^ p{uj) = u{uj) ioi all uj en . (3.4) 

To prove 1)3.11) it is therefore enough to establish that 

hm S (plf^l'^ I z.;(,)) = (3.5) 

for every p G [0, oo). 
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We recall that T^iivi) = J2x=iVx denotes the number of particles. By 
definition ()2.3p one has for a fixed rj & Ql 

™" Z{N,L) • ' ' 

where denotes the product measure W^{r]) = Y[t=i ^iVx)- Thus, using 
()3.3|) . it follows that for every G [0, 0c] 

since n^'^ is absolutely continuous w.r.t. z/^. From ()2.5|) and ()2.6|) we con- 
clude that 

u^{il)Z{<pf = W\'n)(t>'' for -nenL^N and 
{{El = N}) Z{4>f = Z{N, L) 0^ . (3.8) 

To simplify notation we use the shorthand f^^^A) = J^rjeA ^^iv) ^^i a subset 
A of the configuration space and {rj \ ^Liv) = ^} = {^l = N}. Inserting 
in ()3.7p this yields 

S{^,'''^\uj;) = -\oguj;{{J:, = N}). (3.9) 

At this point we use the subadditivity of S, namely if two measures /i, u 
have marginals /ij, i>i, i = 1,2, then 

s{fi\u) > s{^il\ul) + s{fi2W2) . (3.10) 

Therefore for every n G {1, . . . , L} and G [0, 0c] 

S (/^^r-" I 'y;) < log {{El = N}) . (3.11) 

The key point is to maximize z/^ {{^l = N}) by appropriately adjusting 
= 0(p) as defined in ()2.9|) . In the subcritical case, p = N/L < pc, we have 
0(p) < 0c and has exponential moments (see p.6|) ). Then the variance 
cr^ of z/0 is finite and the limit distribution of (S^ — (pL))/{a\/L) is given 
by the normal distribution 7\/'(0, 1) (cf. ()4.1H) ). By the local limit theorem 
(see e.g. [24J) we get in this case z/^^^^ ({^l = [p-^]}) — for large L. 

For p = Pc the decay of is subexponential, since = 0c is the radius of 
convergence of the partition function Z{(j)) in (|2.7p . Thus 

lim ilogz/^,(A;) =0 (3.12) 

and the second moment of z/,^^ could be infinite, leading to a non- normal 
limit distribution (cf. ()4.12|l ). Since the first moment of z/^^ equals pc < oo. 
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by the local limit theorem for non- normal limit distributions (see also j24j ) 
we get the lower bound {{^l = [PcL]}) ^ 1/-^- The supercritical case 
p > Pc where 0(p) = 0c, can be reduced to the critical one via 



^i{{^L = [pL]}) > 

> ^i{{VL = [pL] - [PciL - = [p,(L - 1)]}) = 

= u^^{[pL] - [pe(i^- = [pciL-l)]}) ■ (3.13) 



Both terms decay subexponentially, the first one using ()3.12|) and the second 
one as in the critical case. Thus in all cases we have a subexponential lower 
bound on i/^^^^ ({^l = [P-^]}) and the limit (j3.5p follows for all p G [0, oo) 
from dnUl), 



2m^(/i[^^l'^|z.;(^)) <-2im-^ (3.14) 



To establish ()3.2|) we use the second line of ()3.8|) and immediately get 



for all A^, L. With = (f){N/L) we can use the above estimates, so that the 
second term on the left vanishes in the limit L ^ oo and the assertion ()3.2|) 



Theorem 1 ensures us that the volume fraction of the condensed phase van- 
ishes in the limit L — oo. In principle it could still contain an infinite 
number of sites, and the question remains, how many condensed sites there 
are in a typical configuration. The answer depends on the large-A; behavior of 
the critical distribution ^'^^(fc). From ()3.12|) we know already that it decays 
subexponentially. We will show that for a large class of such subexponential 
distributions the excess particles condense on a single, randomly located site. 

Let C2 be the set of distributions u^^ which have finite second moment 
and for which the integrated tail i^(i,^{{r]i > k}) is heavier than exp[— 
for some a G (0,1/2). Cp denotes the set of distributions with power-law 
tail v^X^) — 1 ^ ^ (2,3], for which the second moment diverges. For a 
detailed description and results on subexponential distributions we refer the 
reader to 0123]. 

Theorem 2. Let z/^^ be in the class C2 U Cp as defined above, with first 
moment pc < 00. Then for the sequence of the corresponding canonical 
measures p^P^^'^ with p > pc one has 



^\ogZ{N,L)-^ 



\og4 ({Si = N}) = logZ(0) - -log0 , (3.15) 



follows. 



□ 




(3.16) 
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The probability that there is a site which contains at least (p — pc)L 
particles converges to one in the thermodynamic limit. But the occupa- 
tion number cannot be substantially larger than (p — pc)L, since almost all 
sites are distributed according to z/^^, as proved in Theorem 1. Therefore 
all excess particles are condensed on a single site in the limit L ^ oo. The 
proof of Theorem 2 uses large deviation results on the asymptotic behavior 
of ({E^ > pL}^ for L ^ oo, which we summarize in the following Lemma 
for our purpose. 



Lemma. Let 'r]i,i]2, ... be i.i.d. random variables with mean zero and 
probability distribution P G C2 U Cp. Then with = ^^^^ one has for 
any p > 



lim 



P^({Si>pL}) 



LP{{r^^>pL}) 



(3.17) 



Proof of the Lemma. In the case P E Cp see |H] (Chapter 1, Corollary 
l.Ll to 1.L3) and for distributions in C2 see jH]- 

The interpretation of the Lemma is that under a distribution in C2 U Cp 
the rare event {T^l > pL} in the limit L ^ 00 is realized by the deviation of a 
single (randomly positioned) site with probability one. Thus for z/^^ G C2UCp 
this is also a typical configuration under the canonical measure p^^^'^, since 
for p > pc the latter is basically given by the grand canonical critical measure 

under the condition {S^ = [pL]}, cf. (j3.19|) . This argument will be made 
precise in the following. 

Proof of Theorem 2. Obviously one has 

pl^^l'^ max r/. > [(p - p,)L]}^ < 1 (3.18) 

for all L E N. To find a lower bound that converges to 1 for L ^ 00 we 
apply the results of the Lemma. We note that for all r] G 0,l,n it is 

max rix>[{p — Pc)L]\) 

"""'({sm* >- - -'^1 }) - 

cf. ()3.6|) . and for L ^ 00 

z/^^^d max [(p- p,)L]}) =L z/^jr^i > [(p- p^L]})(l+o(l)) . (3.20) 



10 



GroBkinsky et al. 



Shifting the expectation value to zero via f]x = Vx — Pc, ^'l •= J2t=i v'x 
using ({Sl = [pL]}) < uj^^ii^L > [pL]}) we get from and 

maxr/.>[(p-p,)L]})> 

Lu^^{{i,>[ip-p,)L]}) 
- -U{^'.>[ip-P.)L]}) ^'-''^ 

By (j3.17|) the righthand side of (|3.2ip converges to 1 in the hmit L ^ oo 
and (|3.16|) is shown. □ 



4 Stationary properties near criticality 



In j7] Evans studies the zero range process with rates 
g,ik) = eik)il + b/k) , 



(4.1) 



where ^(0) = and 6{k) = 1 for A; > 0. He observes that for b > 2 the critical 
density pc < oo. Our goal here is to study the properties of the invariant 
measures for any 6 > 0. The stationary weight for the zero range process 
with rates gt is given by 



1=1 



k\ 



k\T{l + b) 

+ b)k ~ r{i + b + k)' 



(4.2) 



where (a)^ = HiLo^l'^ + denotes the Pochhammer symbol, a E M, k E N. 
The grand canonical partition function of (j2.7|) is 

" (l)fc(l)fc 0^ 
(l + 6)fc k\ ■ 



Z(0)= ,F,{l,l;l + b;<f)):=J2T[ 



(4.3) 



fc=0 



Its radius of convergence is 0c = 1, and 2F1 denotes the hypergeometric 
function [20], which has the expansion 



2Fi{k,k;k + b]. 



T{k + b)T{k-b) 



T{ky 



T{k + b)T{b- k) 



1 + 



r(6)2 f ' l + k-b' 

The particle density ()2.8|) is given by 
2^1(2, 2; 2 + 6; 0) 



\b-k 



1 + 0{1 
+ O(l-0)' 



+ 



(4.4) 



P(0) 



[1 + b) 2Fi(l,l;l + 6;. 



for (h < 1. 



(4.5) 
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In the following we analyze the grand-canonical single site measure of 
()2.6|) in the limit 0/^1, i.e. near the critical density pc- For < oo the 
limit ui is well defined, as discussed in Section 2, and according to Theorem 
1 it is the distribution of the non-condensed phase for supercritical systems 
with N/L = p > Pc- As long as < 1 the distribution has exponential 
moments. For = 1 the exponential tail of disappears and the tail be- 
comes proportional to the weight W{k) as defined in ()4.2|) . Using Stirling's 
formula, the behavior of the weight for large k is given by the power-law 
W{k) ~ r(l + b) k~^. These distributions have moments up to order 6—1. 
Thus different scenarios are encountered as b is varied. 

The case < 6 < 1: 

For b < 1 the leading order in the asymptotic expansion for Z and p is given 

by 

z(0)~r(i + 6)r(i-6)(i-0)''-i ^ oo , 

- (1 + _ 5) (1 - 0)-^ - Pc = oo , (4.6) 

as ^1. For every density, the stationary distribution in the limit L ^ oo 
is given by the grand-canonical measure (see Figure^a)). The probability 
to have a fixed number of particles on a given site vanishes with increasing 
density, as is shown for the example of an empty site in Figure ^b). Thus in 
the limit there is an infinite number of particles on every site with probability 
one, as it should be for homogeneous systems with p —>■ oo. 

For 6 = 1 this picture does not change qualitatively, except for the loga- 
rithmic corrections 

log(l - 0) 

^(0) = oo , 



^^'^^^ (0-1)4(1-0) -^ ^ ^^=^' ('-^^ 

as /" 1. 



The case 1 < 6 < 2: 

For 1 < 6 < 2 the leading order terms change, and 

z(0)^^ + r(i + 6)r(i-6)(i-0ri - zi^) = ^^ 

p(0) ^ 0(6 - l)r(6)r(2 - 6) (1 - <pf~' ^ = oo , (4.8) 

as y I. As before, for 6 = 2 the first order terms have logarithmic correc- 
tions but the qualitative behavior does not change. In particular, pc = oo 



12 



GroBkinsky et al. 




12 14 



Figure 1: (a) Fugacity (/> as a function of the particle density for several values 
of b. Equivalently, current j as a function of p up to the scale factor m{p) (see 
(|2.11|) ). (b) Probability of an empty site, j^0(p)(O), as a function of the density p 
for several values of b. 



and the stationary distribution is described by the grand-canonical ensemble 
for every density p (see Figure Ufa)). 

However, somewhat surprisingly, the character of this distribution for 
large p differs from the case b < 1. Since Z{1) < oo, ui is well defined and 
there is a non-zero probability to have a fixed number of particles at a given 
site, 

1 _ b-1 

^^^^ ^ r{b) [b - 1) k-^ for large k. (4.9) 
z(l) 

For example the probability of an empty site, given by z^0(p)(O) = 1/Z(0(p)), 
decreases monotonically with 0, i.e. with increasing density p. In contrast to 
the case b < 1, it does not vanish in the limit p — > cxd, however, it reaches the 
non-zero value z^i(O) = {b — l)/b (see Figure ^b)). So no matter how large 
the density, the fraction of empty sites in a typical configuration is always 
greater than (6 — l)/b. 

Distributions with power-law tails are well studied (see e.g. j27j and ref- 
erences therein). A typical configuration for this stationary distribution, i.e. 
a set of L i.i.d. random variables 7]^ drawn from z/i, is known to have a hi- 
erarchical structure. The n-th largest value of the set {rji, . . . ,riL} scales as 
(r(6 — l)L/nY^'^''~^\ which holds for every 6 > 1. In our particular case 
1 < 6 < 2, this means that the particle number T.^ also scales as L^/^^^^^ and 
thus grows faster than the number of summands L. Therefore the particle 
density T^l/L diverges as L^^"''^/*^''"^^ and the highest occupied site contains a 
nonzero fraction of the particles in the system. This hierarchical structure of 
typical configurations can be understood as a precursor for the condensation 
phenomenon to be discussed in the next part. 



^i(0) = 
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The case b > 2: 

In this case b is large enough so that besides the normahzation also the first 
moment of the grand canonical distribution converges in the limit /" 1: 

Z(0) ^ - — -4— ^ (1 - 0) - Z{1) 



6-1 (6- 1)(6-2) ' ^ ' 6-1 ' 

p{<i^)^^+<Pib-mm2-b){l-<j^f-' - Pc = ^, (4.10) 

as y 1. We note that for 6 > 3 also the second moment of the 
distribution ui exists and the number of particles satisfies the usual central 
limit theorem 

hm z/f fei < < 6) = r G{0 dt (4.11) 

where G denotes the Gaussian probability density with zero mean and unit 
variance. The density for 6 > 3 is of order p{(f)) = 1/(6 — 2) + 0(1 — 0) and its 
first derivative is finite at = 1 and given by p'(l) = (6 — 1)^/((6— 3)^(6— 2)^) 
(see Figure m^a)). 

As explained already in Section 3, the most occupied site contains of 
order L^/^'^~^') particles, and for 6 < 3 this fluctuation is larger than a/L. 
Therefore the scaling limit leads to a self-similar distribution, which is given 
by the completely asymmetric Levy distribution i^(b-i) (for details see 

or m), 



lim 

L-^oo 



( y , \ 

^1 ^ ^ N 1/(6-1) ^ ^2 

\^ (6r(6-l)L) J 



\(,_i)(Orfe (4.12) 



With ()4.10j) we have p'(l) = oo for 6 < 3, leading to a differentiable function 
0(p), as shown in Figure [H^a). 



5 Dynamics of the condensation 

The stationary distribution investigated so far carries no information on the 
kinetics of the condensation. A natural set-up is to start with particles 
uniformly distributed at the supercritical density p > pc- In the beginning 
the excess particles condense at a few random sites. Such a site containing 
many excess particles is called a cluster site. Thus there are several clusters 
which are essentially immobile as will be discussed below. On the remaining 
sites, called bulk sites, the distribution relaxes to ui. With increasing time 
the larger clusters will gain particles at the expense of the smaller ones. 
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causing some of the clusters to disappear. Eventually only a single cluster 
containing all excess particles survives, which is typical for the stationary 
distribution, as has been discussed already in Section 3. 

In the following we will study the kinetics of condensation and its depen- 
dence on system parameters in detail. Compared to Section 4 the jump prob- 
abilities play an important role for the dynamical properties of the system. 
We will focus on nearest-neighbor jumps which are either totally asymmet- 
ric, i.e. particles only jump to the right with Pa{y) = 6{l,y), or symmetric, 
i.e. Ps{y) = 1,1/) + 6{l,y)) /2. In contrast to previous sections we are 
not able to rigorously prove our statements, but use heuristic considerations 
which are corroborated through comparison with simulation data. 

5.1 Cluster formation 

To distinguish cluster and bulk sites, we define a site x G {1, . . . , L} to be 
a cluster if it contains a macroscopic fraction of the excess particles rj^ > 
a{p — pc)L. The prefactor a G (0, 1] is rather arbitrary but for simulations it 
is important that the clusters are well separated from the bulk fluctuations, 
which are of order {T{b — l)Ly^^^~^\ Since fluctuations grow only sublinearly 
with L this separation is clearly guaranteed in the limit L —>■ oo, and for 
finite systems it holds for sufficiently large L depending on a and p. In our 
simulations we choose a = 1/40, requiring system sizes of about 200 sites 
minimum for the values of p and b considered. Let n{t) be the number of 
cluster sites at time t and mj(t), i = 1, . . . , n(t) be the size of the i-th largest 
cluster, i.e. mi(t) > ... > mn(t)(t). These quantities depend also on the 
system parameters b, L, and p. 

By definition a typical cluster has a size of order a{p — Pc)L and so there 
are of the order of 1/a cluster sites. The time scale for the formation of 
such clusters is very roughly estimated as follows: 0{L) particles have to 
move a distance of order L to form the cluster. So in the asymmetric case 
the time scale for cluster formation is 0{L'^). The dependence on p and b 
is not so obvious, since the bulk has not yet relaxed to vi and the speed of 
particles still changes. In the symmetric case the time scale is 0{L^), since 
the particles diffuse without a drift. 

In Figure|2fa), in the totally asymmetric case, we plot the average number 
of particles in the condensed phase, normalized by the number of excess 
particles. 



fit) 



{p - Pc)L 



(5.1) 



The time axis is scaled proportional to Ta = {p — pc)'^L'^/b (cf. Section 5.2), 
since this choice gives the best data collapse when p and b are varied. On 
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t / Ta t / Ta 

Figure 2: (a) Average fraction (/(t)) of particles in the condensed phase, (b) 
average number of cluster sites {n{t)). Both are plotted as a function of time in 
units Ta = {p — PcYlP' /h (cf. (|5.6() ') for different values of p, h and L in the totally 
asymmetric case. Symbols: for 6 = 4, p = 5, L = 320(0), 640(A), 1280(n), 
2560(x), for 6 = 4, p = 3, L = 1280 (*), for 6 = 5, p = 5, L = 2560 (■). 



this time scale most excess particles become trapped in a cluster and the 
bulk relaxes to z/i. In Figure Efb) the average number of clusters is plotted 
as a function of time. The number of clusters grows for a short time and 
then starts decreasing again. 

5.2 Coarsening 

Once the bulk has relaxed to z/i, each bulk site loses particles at the average 
rate {gb)vi = 1- In the asymmetric case this results in a particle current 
i = (p = 1 (see ()2.11|) and Figure H]). On top of that, excess particles are 
exchanged between clusters. The bulk can be seen as a homogeneous medium 
where the excess particles move, and the cluster sites as boundaries where 
they enter and exit. A cluster of size m > loses excess particles with rate 
5f(m) — 1 = h/m and gains particles from neighboring clusters. Since this 
rate decreases with increasing cluster size, smaller clusters lose particles to 
the larger ones. 

To quantitatively describe this coarsening process we study the normal- 
ized mean cluster size fh{t) = f(t)/n{t) as a function of time for large system 
sizes L. The ensemble average of this quantity (denoted by (...)) is expected 
to grow according to a scaling law 

(m(t)) ~ , (5.2) 

with a scaling exponent f3 To estimate this exponent we notice that the 
time scale for the coarsening process is determined by two factors: Firstly, 
the rate at which excess particles leave a cluster of size m and enter the bulk. 
Secondly, the typical time of such a particle to reach the neighboring cluster. 
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The mobility of excess particles in the bulk is characterized by the average 
exit rate from an occupied site, 

^ z/i(fc 7^ 0) 1 - z/i(0) V b J 

as follows from (j2.11|) and (|4.9|) . The excess particles perform a random 
walk in the bulk, which is either biased or unbiased, depending on the first 
moment of the jump probabilities. The time it takes to reach the next cluster 
at a distance of order 0{L/n) = 0{m/{p — pc)) is given by the mean first 
passage time for a random walk [21]. In the following we distinguish between 
asymmetric and symmetric jump probabilities, to estimate the time scale for 
a particle to leave a cluster and the one for reaching the neighboring cluster. 

In the totally asymmetric case excess particles leave a cluster of size 
m with rate b/m and move towards the right neighboring cluster without 
returning, so the time to lose one particle scales like 0{m/b). Since the mean 
first passage time of a biased random walk is proportional to the distance, 
this particle spends a time 0(m/(6(p — Pc))) in the bulk, because b is the 
speed of the particle. Thus the typical times for exiting a cluster and entering 
the next neighbor are of the same order in L and there are 0{l/{p — pc)) 
excess particles in the bulk. The coarsening time scale, determined by the 
typical time for a cluster to lose all m particles, is thus proportional to 

2 

771 

ta{m) := — , i.e. Pa = 1/2 . (5.4) 

So for the asymmetric case the mean cluster size is predicted to grow like 
(m(t)) ~ (6t)V2. 

In the case of symmetric jump probabilities ps, excess particles perform an 
unbiased random walk in the bulk with diffusion constant b. Thus the mean 
first passage time to reach the next cluster is proportional to the square 
of the distance, i.e. 0{rn? /{b{p — PcY))- In contrast to the asymmetric 
case, it is very likely that particles return to the cluster they left. The 
probability that they do not return but reach the neighboring cluster, which 
is the relevant event for coarsening, is inverse proportional to the diffusion 
distance, i.e. C((p — Pc)/Tn). So the typical time of a particle to leave a 
cluster is 0{rn? /{b{p — Pc))) and as before there are 0(l/(p — pc)) excess 
particles in the bulk. The coarsening time scale is therefore proportional to 

Tfi 771 

ts{77i) := ta{77i) = — , i.e. = 1/3 . (5.5) 

(P - Pc) b[p - Pc) 

Thus the mean cluster size is predicted to grow like {7h(t)) ~ ((p — pc)bty^^. 
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Figure 3: Double-logarithmic plot of mean cluster size as a function of time. 
Data collapse is achieved by using the appropriate time scales r^, given in (|5.6j) . 
and dividing the cluster size by the number of excess particles {p — pc)L. The 
straight line indicates the predicted slope, (a) /? = 1/2 in the asymmetric case, 
(b) /? = 1/3 in the symmetric case. Symbols: (a) 6 = 4, p = 5 and L = 320(0), 
640(A), 1280(n); (b) 5 = 4, p = 3 and L = 160(0), 320(A), 640(n). 



In general, the time scale for the coarsening regime is determined by the 
largest clusters with size C((p — Pc)L) and is thus of order 

Ta:=ta{ip- Pc)L) = {p-p,fL^/h 

Ts:=ts[{p-Pc)L) = (p-p^yiyb (5.6) 

for asymmetric resp. symmetric jump probabilities. The growth exponents 
for the two cases are confirmed by simulations and shown in Figure Efa) 
for the totally asymmetric and Efb) for the symmetric jump probabilities. 
Using the time scale resp. and normalizing (m) by the number of excess 
particles {p — pc)L the data for different system sizes collapse. The measured 
growth exponents from these data are 

Pa = 0.514 ± 0.005 and /5, = 0.334 ± 0.004 , (5.7) 

which agree with the above predictions. Independently from us, these expo- 
nents have been obtained in |T|]j by numerical simulations. 

Note that the clusters coarsen on the same time scale as they nucleate. 
However, looking at the number of time steps in Figures El and El the time 
it takes to nucleate clusters is by a factor of 10 shorter than the coarsening 
regime. 

5.3 Saturation 

Eventually all clusters except for two will have disappeared and finite size 
effects become dominant. The scaling law ()5.2j) is no longer valid in this 
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Figure 4: Average size of the three largest clusters (mj(t)), z = 1,2,3 as a 
function of time in the totally asymmetric case (a) and the symmetric case (b). 
Data collapse is achieved by using the appropriate time scales , given in (|5.6|) , 
and dividing the cluster size by the number of excess particles {p — pc)L. Symbols: 
(a) 6 = 4, /9 = 20 and L = 80(0), 160(A), /o = 40 and L = 80(n), 160(x). (b) 
6 = 4, p = 20 and L = 40(0), 80(A), p = 40 and L = 40(n), 80(x). 



regime, since the mean cluster size fh saturates towards its limiting value. 
The two clusters exchange particles until one of them vanishes and the system 
has reached its stationary state where all excess particles are concentrated at 
a single cluster site. In Figure ID we plot the average size of the three largest 
clusters (mi(t)), i = 1, 2, 3 normalized by {p — pc)L. Note that the coarsening 
regime ends at latest when the third largest cluster has disappeared, and thus 
takes only about a tenth of the total equilibration time. In the following, 
we focus on the totally asymmetric jump probabilities, but the symmetric 
choice would lead to an effective evolution equation of the same form. 

Let M = mi + m2 be the total number of particles at the two largest 
cluster sites. On the time scale t* = t / {{p—pc)L/b) the two clusters exchange 
single particles with effective rates (p — pc)L/mi, i = 1,2 (see discussion in 
Section 5.2). The fluctuations of M on this time scale are only 0{1). Thus it 
is M = {p — pc)L + 0{l), since the bulk is relaxed to ui and all other clusters 
have disappeared. Let g(m, t*) be the probability of having m = 0, . . . , M 
particles at one cluster site and M — m at the other one. The dynamics is 
then governed by the effective master equation 



d 

—q{m,t*) = -q{m,t* 



e{m) _^ e{M -m) 



m/M 1 - m/M _ 



+ 



0{m) , e(M-m) , , 

^('"-M-) ^_(„'_i)^,, +.(.n.+ l,r),^A^. (5.8) 

The exchange rates on the right hand side only depend on the rescaled vari- 
able m/M and not on the system parameters p, L and h. For large M, m/M 
varies on the time scale t* / M 'r^t / {{p — pc)'^L'^/b), confirming that fl5.()|l 
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is the appropriate time scale for the two-cluster situation. Therefore the 
plots in Figure m are independent of the system parameters. However, in the 
following discussion we stick to the time scale t* and the discrete variable m. 

For any initial condition the solution of ()5.8|) tends to the inverse binomial 
distribution q*{m) ~ 1 / (^) . It is symmetric around m = M/2 and for small 
m we have q*{m) = 0{M~'^). Thus the two extreme occupation numbers 
m = resp. m = M are the most probable ones and in the limit L — oo both 
have probability 1/2, consistent with the results of Section 3. For m = aM, 

a e (0, 1) it is g*(m) = C^\/M(^a"(l - a)^"")*^^ using Stirling's formula. 

Thus, in the stationary state, the typical time for a macroscopic fluctuation 
of the cluster size diverges exponentially with the system size L. 

To study the relaxation dynamics, we write ()5.8p in the canonical form, 
using the discrete derivative Vm f{rn) := /(m + 1) — f{m), 

— g(m, t*) = -Vm (a{m) q{m, t*)^ + (^rf(m) g(m, t*) 



aim) 



1 M 2m/M - 1 



(1-m/M) m m/M{l-m/M) 



1/1 M\ _ 1 

~ 2 V(l - m/M) ^m)~ 2m/M(l - m/M) ' ^ ' ' 



For ease of notation, we take < m < M and ignore the boundary terms. 
Note that ()5.9|) is symmetric around m = M/2. It describes diffusive motion 
in a double well potential with drift a(m) and diffusion coefficient d{m), with 
the slightly unusual feature that the minima of the potential are located close 
to the boundaries at m = 1 and m = M — 1. 

The master equation must be supplied with a suitable initial condition 
g(m, 0), which, since resulting from a complex coarsening process, is not 
readily available. A crude estimate can be found by noting that g(m, 0) 
is roughly proportional to the lifetime of the occupation number m. It is 
determined by the inverse exit rate for m taken from Equation (|5.8p . and 
thus we expect g(m, 0) ~ 6m/M (1— m/M). This is a symmetric single hump 
distribution with mean M/2 and standard deviation M/{2\/5) ~ 0.22M. 
Comparing with the simulation data at the time when the third largest 
cluster has just disappeared, m^{t)/M < 0.01, we indeed find a single hump 
distribution with mean M/2 and standard deviation 0.166M. When solving 
Equation (j5.9j) with this initial distribution, the expectation of the larger 
cluster size, given by 1/2 + (|m/M — l/2|)g(m,t), is indistinguishable from 
(mi (if:)) in Figure El Note that, except for the time scales. Figure Ufa) and 
I3fb) are almost identical, confirming that the effective master equation for 
the symmetric case is of the same form as ()5.8j) . 
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